Prolate dominance of nuclear shape caused by a strong interference between the 
effects of spin-orbit and P terms of the Nilsson potential 
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The origin of the dominance of prolate shapes over oblate ones of the ground states of atomic 
nuclei is investigated with the Nilsson-Strutinsky method. The number of prolate nuclei among all 
the deformed even-even nuclei is calculated as a function of the strengths of the spin-orbit and the 
terms of the Nilsson potential. The latter simulates a square-well like radial profile of the mean 
potential. The ratio of prolate nuclei is 86% with the standard strengths corresponding to the actual 
atomic nuclei. By weakening the spin-orbit potential, the ratio oscillates strongly, having a local 
minimum value of 45% with a half of the standard strength and a local maximum value of 78% 
without the spin-orbit potential. 
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A basic question in nuclear physics is why atomic nuclei have a strong tendency to deform into prolate shapes than 
oblate ones. 

^5 , Since the early days of the discovery of nuclear deformation it has been usually believed that the nuclear 

■ deformation can be ascribed to the shell structure of nucleon's single-particle spectrum. One might suspect the 
, existence of some unknown simple and direct correspondence between the prolate dominance and a feature of the 

hamiltonian, e.g., a specific term of the elementary nucleon-nucleon interaction, in an analogous fashion as the tensor 
" ^ , force causes a mixture of d wave in the wavefunction of a deuteron. However, it is probably sufficient at the present 
^ I ' stage to confine the scope of the investigation to the mean single-particle potential, through which most of the possible 
\^ ', causes affect the deformation. 

• There are two causes which favor prolate shapes not through the shell effect. One is the Coulomb repulsion between 
' protons, which tends to deform the nucleus into an elongated shape rather than a flattened shape. This effect is, 
' however, important only in heavy nuclei while the prolate dominance is present already in middle-weight nuclei. The 
other, argued by Zickcndraht 1^, is the difference of the volume element of the collective coordinates between prolate 

■ and oblate shapes, which can be identified with the difference of the available configuration space in spherical shell 
''pi \ model calculations. In mean-field approaches, this effect corresponds to that of an angular momentum projection into 

' I ■ spin zero states. However, it does not seem to be essential because the prolate dominance can be reproduced without 
O , the projection in shell-correction ^ and mean-field methods. 

2 ' Let us mention three kinds of potentials which give rise to shell effects favoring prolate shapes. 

' The first one is the anisotropic harmonic oscillator, which is the most simple approximation used for the nuclear 
mean-field potential. Concerning the sd shell nuclei, Bohr and Mottelson stated that prolate (oblate) shape is preferred 
k>( in the beginning (end) of the major shell filling due to the strong shape driving effect of the particles (holes) in the 
n = I orbital 1). This seems to suggest an equal number of prolate and oblate nuclei. According to Ref. it is 

\ the spin-orbit potential which breaks this even situation by weakening the oblate-shape shell effect. A more general 
argument given by Castel et al. |^ is that the summation of the single-particle energies of an isotropic harmonic 
oscillator is decreased by extending one axis and shrinking the other two axes under volume conservation condition 
neglecting detailed effects of the Pauli principle. (Their argument seems to apply only to the harmonic oscillator 
potential contrary to their insist.) Therefore, harmonic oscillator potentials are expected to favor prolate shapes. A 
quantitative estimation of this effect is one of the aims of our study. 

The second kind of potentials are those with square-well like radial profile. The nuclear mean potential resembles 
the Woods-Saxon potential [Q, which is in between a square well and a harmonic oscillator. Frisk found that such 
radial dependence is an origin of the prolate dominance from an analysis of classical periodic orbitals in an ellipsoidal 
cavity. By considering the volume conservation, he showed that the shell effect at the Fermi surface changes strongly 
in the prolate side while stays almost constant in the oblate side as a function of the magnitude of deformation. 
Consequently, if the spherical shape is unstable, oblate shapes are equally unstable but there must be a more stable 
state in the prolate side. 

The third one is the spin-orbit potential, which is indispensable for the reproduction of the spherical magic numbers 
and is an important component of the nuclear mean potential. Its relation to the prolate dominance is suggested 
from an extensive Skyrme-Hartree-Fock calculation The energy difference between prolate and oblate minima 
exhibits a clear and abrupt change of behavior between Z,N < 40 and Z,N > 50 where Z and are the numbers 
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of protons and neutrons, respectively. In the former region prolate and oblate solutions appear evenly in the ground 
state, while in the latter region the oblate solutions have systematically higher energies than prolate ones. Between 
the two regions, the character of major shells change from the harmonic oscillator type to the Mayer- Jenzen type, the 
latter of which includes a high-j intruder in each major shell due to the spin-orbit potential. This parallelism suggests 
that the spin-orbit potential plays an essential role in giving rise to the prolate dominance. 

In this paper we will examine how the situation of the prolate dominance changes when the radial profile of the 
potential and the strength of the spin-orbit potential are different from those of actual nuclei. For this purpose, we 
employ the Nilsson-Strutinsky method which is a convenient and well-established method to reproduce nuclear 
shapes. The single-particle potential of the method is called the Nilsson or the modified oscillator potential and is 
expressed as, 



The first term stands for an anisotropic harmonic oscillator potential, where the frequencies lo±_ and ljJ\\ are expressed 
as functions of a quadrupole deformation parameter £2, 



deformation potential. The third term is a spin-orbit potential, in which orbital and spin angular momenta are 
expressed as 1 and s, respectively. Subscript t means the usage of the stretched coordinates. The fourth term is called 
the P term or potential hereafter. The Woods-Saxon type radial dependence of the potential is approximated by this 
term including the square of the orbital angular momentum. The standard values given in Table 1 of Ref. ]lO[ | are 
used for the parameters km and /iAr, which are dependent on the total of the oscillator quanta, N . The factors fu 
and fu are introduced in this paper to modify the standard potential, which is restored by putting fis — fii — 1- A 
convenient feature of the Nilsson potential for our study is that the spin dependent and independent potentials can 
be changed independently unlike in the Woods-Saxon potential or the relativistic mean- field model pT[ |. 

We have utilized a program provided by Y. R. Shimizu ||T^. It is based on the NIGRA code but is simplified 
for non-rotating axially symmetric states, which makes the calculation much faster. The pairing correlation is active 
for single-particle levels within ±I.2?iwo from the Fermi level, while the strengths of the pairing force are determined 
such that the pairing gap for smoothed level density becomes A = 12, A~^/^ MeV. The parameters of the macroscopic 
part are = 17.9439 MeV, Ks = 1.7826, and i?c = l.224QA^/^ fm. See Ref. ||l|l for the details of the model. 

Calculations with the above model have been done in the following way. We choose the values of reduction 
factors /;/ and fis. For each combination {fii,fis), we calculate the potential energy surface (PES) curve versus £2 
(—0.5 < €2 < 0.5, with £4 optimized in —0.16 < £4 < 0.16 for each £2) for all the even-even nuclei with 8 < Z < 126 
and 8 < < 184 and between proton and neutron drip lines predicted by the Bethe-Weizsacker mass formula 0. The 
number of nuclei thus included is 1843. We neglect the possibility of triaxial deformations since non-axial shapes are 
very rare for even-even nuclei 0|. The reduction factors are taken from a square area — 1 < < 1 and —1 1. 
We have taken the sampling spacing to be Afu = A/;^ = 0.125, which has resulted in a computation time of a few 
months with a latest personal computer. 

For each PES curve, we have to label the shape of the ground state as prolate or oblate. One has to be very careful 
in generating an algorithm for this purpose. Well-deformed nuclei usually have both prolate and oblate minima and 
thus each minima can be labeled without ambiguity. On the other hand, transitional nuclei often have several shallow 
minima in a large valley extending from oblate side to the prolate side. In such a situation, it is not meaningful 
to discuss which minima has the lowest energy. After examining a large number of PES curves, we have decided 
to consider only those nuclei which have both well-developed prolate and oblate minima. The practical procedures 
we finally adopted are as follows: (1) Draw a smeared PES curve obtained through a convolution with a weight 
fimction exp[— (A£2/0.05)^]. (2) Separate the original (i.e., before the smearing) curve into valleys by regarding local 
maxima of the smeared curve as the "water-shed" . (3) Taking up the minimum of the original curve in each valley, if 
£2 < — 0.05(> 0.05) at the minimum and £2 < 0.1(> —0.1) at the right (left) end of the valley, regard the minimum as 
representing an oblate (prolate) solution. (4) If a nucleus has both oblate and prolate solutions satisfying the above 
criteria and the deeper one is the oblate (prolate) one, count the nucleus as an oblate (prolate) nucleus. Denoting thus 
counted number of oblate (prolate) nuclei with No {Np), we define the ratio of prolate nuclei as Rp = Np/{Np + No). 
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Rp may take values from to 1. The denominator Np + No is about 900 on the average. Note that the smeared curve 
is used only to divide the curve into valleys and it does not affect the energy or the deformation of the minima. 

Fig. H shows the ratio of prolate nuclei Rp as a function of the reduction factors {fiufis) by means of contours for 
Rp and symbols for the locations of local maxima (triangles) and minima (squares). Fourth order polynomials in 
or fis are used for the interpolations to draw the contours and locate the extrema. 

The standard nuclear potential corresponds to top-right corner of the figure, where Rp takes the largest value in 
the entire square area. The value is Rp = 86%: Among about 900 even-even nuclei having both prolate and oblate 
minima, 86% are prolate in the ground state. One can say that the prolate dominance is reproduced with the standard 
Nilsson potential. 

Our result is also in qualitative agreement with a calculation for metallic clusters p^ ], in which prolate ground 
states are found to be roughly twice as many as oblate ones in the framework of the jellium model with infinite square 
well potential. The corresponding region in Fig. is //, = and where Rp -70 - 80%. 

The minimum value of Rp is obtained at (/;;, //s) = (—1, —0.125), where Rp — 40%, i.e., 60% of the deformed nuclei 
are oblate. The increasing trend of Rp as a function of fu along //s = line implies that the attractive (repulsive) 
term favors prolate (oblate) shapes. This result confirms the view of Frisk. 

On the other hand, the spin-orbit term cannot be regarded as favoring either prolate or oblate shapes, because Rp 
behaves roughly symmetrically between positive and negative values of . The most conspicuous fact concerning the 
spin-orbit term found in our study is a very strong interference with the P term. In Fig. |l|, by moving down from the 
top-right corner along a line = 1, Rp takes on 86, 45, 78, 44, and 81% for fis=^, \, 0, — i, and —1, respectively. 
One can see two periods of oscillations in — 1 < /js < 1. Weakening the spin-orbit term by 50% moves the ratio Rp 
from the highest peak at fis = 1 to the bottom of a deep valley at fis = ^, where there are more oblate nuclei than 
prolate ones. A complete disappearance of the spin-orbit term moves the ratio to another high peak at fis = and 
recovers the prolate dominance. Combination of the two terms produces a situation which is beyond expectation from 
the independent effects of each term. 

A prolate dominance as high as 80% is realized only for restricted combinations of the strengths of the two terms. 
It may not be a mere coincidence that the potential of actual nuclei matches one of such rare combinations. The 
same kind of subtle balance between the two terms has been discussed concerning the pseudo-spin symmetry [p^[-p3[ , 
which holds when /i^r of the Nilsson potential is | while the standard values of are between 0.5 and 0.6. 

The prolate dominance can be related to the pseudo-spin symmetry by extending Frisk's view to particles with 
spin: An attractive term can cause prolate dominance if the spin is decoupled from the orbital motion. The prolate 
dominance occurs at = ±1, and fu = 1. The real spin is decoupled at fis = 1 while the pseudo-spin is decoupled 
at fis = 0. The point at (/;;, fia) = (1, —1) might correspond to a similar situation in which another kind of spin-like 
quantity is decoupled. There are opinions that the pseudo-spin symmetry has a relativistic origin. If they would be 
approved finally, the prolate dominance of nuclear deformation might also be related with some relativistic aspect of 
the atomic nuclei. 

Let us mention some other results from our calculations before concluding the paper. (1) The inclusion of optimized 
hexadecapole deformation has a tendency to favor prolate shapes. By setting £4 = 0, Rp is reduced from 86% to 82% 
at {fu, /;s) = (l,l). An average of Rp over all the combinations of the reduction factors is decreased from 60% to 55%. 
(2) An almost pure harmonic oscillator potential (i.e., with fis = fu = ~ ^ and weakened pairing) produces Rp = 
55%. This is an quantitative estimation of the tendency of prolate preference predicted by Castel et al. 

In summary, a strong interference is found between the effects of the spin-orbit and the l^ terms of the Nilsson 
potential. The ratio of prolate nuclei among well-deformed even-even nuclei is more than 80% by using the standard 
strengths for the two terms. Multiplication of ±1 or to the strength of the spin-orbit term does not change the 
situation of prolate dominance. On the other hand, when the strength is multiplied by ±i, the ratio is less than 50%, 
i.e., there are more number of oblate nuclei than prolate ones. The emergence of prolate dominance for restricted 
combinations of the strengths of the two terms is in a parallel situation with the decoupling of real or pseudo spins 
from the orbital motion and can be understood by extending Frisk's view. 

We are planing to study possible changes due to (1) reductions of pairing force strengths and (2) a replacement 
of the Nilsson potential with the Woods-Saxon potential. It is also an interesting question to which region of Fig. |l| 
neutron-rich unstable nuclei, which are waiting for experimental studies, correspond. In such nuclei, both terms are 
expected to be more or less weakened compared with those for stable nuclei |Q . 
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FIG. 1. The ratio of prolate nuclei Rp. The abscissa and the ordinate are the reduction factor of the strength of the 
potential (/;;) and that of the spin-orbit potential (fis) relative to the standard values. Contours are for Rp = 45, 50, 55, • • • 80%. 
Thick curves are for i?p=50%. Solid triangles (squares) indicate the location of local maxima (minima). 
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